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We investigate the critical properties of the Ising 5=1/2 and 5 = 1 model on (3,4,6,4) and 
(3 4 , 6) Archimedean lattices. The system is studied through the extensive Monte Carlo simulations. 
We calculate the critical temperature as well as the critical point exponents y/v, /3/v and v basing on 
finite size scaling analysis. The calculated values of the critical temperature for 5 = 1 are ksTc/J = 
1.590(3) and k B T c / J = 2.100(4) for (3,4,6,4) and (3 4 , 6) Archimedean lattices, respectively. The 
critical exponents fi/v, y/v and 1/v for 5 = 1 are p/v = 0.180(20), y/v = 1.46(8) and 1/v = 0.83(5) 
for (3,4,6,4) and 0.103(8), 1.44(8) and 0.94(5) for (3 4 ,6) Archimedean lattices. Obtained results 
differ from the Ising 5=1/2 model on (3, 4, 6, 4), (3 4 , 6) and square lattice. The evaluated effective 
dimensionality of the system for 5 = 1 are D eB = 1.82(4) for (3, 4, 6, 4) and D bB = 1.64(5) for (3 4 , 6). 
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I. INTRODUCTION 

The Ising model [l], Q remains probably the most 
cited model in statistical physics. Today, the 75/ Web 
of Knowledge abstracting and indexing service returns 
over eleven thousands records for the query on "Ising" 
for time span from 1996 to 2010. For Ins-pec database 
(for years 1969-2010) this number is almost doubled and 
reaches 19 thousands for Scopus database (for data range 
1960-2010). The latter means that during the last half 
of century ks 380 papers refer to the Ising model every 
year. The Google search engine indicates over 279 thou- 
sands web pages which contain "Ising model" phrase. 

The beauty and the popularity of this model lies in 
both its simplicity and possible applications from pure 
and applied physics, via life sciences to social sciences. In 
the way similar to the percolation phenomenon, the Ising 
model is one of the most convenient way of numerical 
investigations of second order phase transitions. 

In the simplest case, the Ising model may be used to 
simulate the system of interacting spins which are placed 
at the nodes of graphs or regular lattices. In its basic ver- 
sion only two values of the spin variable are available, i.e. 
S = — \ and S = +|. This is the classical Ising S = | 
model. For a square lattice this model defines the univer- 
sality class of phase transitions with analytically known 
critical exponents which describe the system behaviour 
near the critical point. The critical point separates two 
— ordered and disordered — phases. 

One of possible generalisation of the Ising model is to 
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enlarge the set of possible spin values (like in the Potts 
model 0,0]). The Ising 5=1 model corresponds to three 
possible spin values, i.e. S G { — 1,0, +1}, Ising S = § 
allows for four spin variables S G {±|,±-|}, etc. The 
Ising S ^ \ model on various networks and lattices may 
form universality classes other than the classical square 
lattice Ising model. 



The spin models for 5=1 were extensively studied by 
several approximate techniques in two and three dimen- 
sions and their phase diagrams are well known (5l4*l~ll|. 
The case S > 1 has also been investigated according to 
several procedures [l2|-[l8j. The Ising model 5 = 1 on 
directed Barabasi-Albert network was studied by Lima 
in 2006 [l9| . It was shown, that the system exhibits first- 
order phase transition. The result is qualitatively differ- 
ent from the results for this model on a square lattice, 
where a second-order phase transition is observed. 



In this paper we study the Ising 5 = 1 model on 
two Archimedean lattices (AL), namely on (3, 4, 6, 4) and 
(3 4 ,6). The topologies of (3,4,6,4) and (3 4 ,6) AL are 
presented in Fig. [TJ Critical properties of these lat- 
tices were investigated in terms of site percolation in Rcf. 
[20j . Topologies of all eleven existing AL are given there 
as well. Also the critical temperatures for Ising 5 = | 
model [2l| and voter model [22| on those AL were esti- 
mated numerically. 



Here, with extensive Monte Carlo simulations we show 
that the Ising 5 = 1 model on (3,4,6,4) and (3 4 ,6) 
AL exhibits a second-order phase transition with criti- 
cal exponents that do not fall into universality class of 
the square lattice Ising S = h model. 



2 




XXXXKXXXX 



XXXXXXX 







FIG. 1: Topology of (3,4,6,4) [left] and (3 4 ,6) [right] AL. 



II. MODEL AND SIMULATION 

We consider the two-dimensional Ising S = 1 model 
on (3, 4, 6, 4) and (3 4 , 6) AL lattices. The Hamiltonian of 
the system can be written as 



N N 



(i) 



2—1 j>i 



where spin variable Si takes values —1, 0, +1 and deco- 
rates every N = 6L 2 vertex of the AL. In Eq. {1} J is 
the magnetic exchange coupling parameter. 

The simulations have been performed for different lat- 
tice sizes L — 8, 16, 32, 64 and 128. For each system 
with N = 6L 2 spins and given temperature T we per- 
formed Monte Carlo simulation in order to evaluate the 
system magnetisation m. The simulations start with a 
uniform configuration of spins (Si = +1, but the results 
are independent on the initial configuration). It takes 
10 5 Monte Carlo steps (MCS) per spin for reaching the 
steady state, and then the time average over the next 10 5 
MCS are estimated. One MCS is accomplished when all 
N spins are investigated whether they should flip or not. 
We carried out N Tnn = 20 to 50 independent simulations 
for each lattice and for given set of parameters (N,T). 
We have employed the heat bath algorithm for the spin 
dynamic. 

We evaluate the average magnetisation M, the suscep- 
tibility and the magnetic 4-th order cumulant U : 



k B T 



M(T,L) = (\m\), 



X (T,L)=N((m 2 )-(\m\) 2 ), 



U(T,L) = 1 



(m A ) 
3(H) 5 



(2a) 



(2b) 



(2c) 



where m = Si/N and ks is the Boltzmann constant. 
In the above equations (...) stands for thermodynamic 
average. 

In the infinite- volume limit these quantities ([5]) exhibit 
singularities at the transition point Tc- In finite systems 



the singularities are smeared out and scale in the critical 
region according to 



M = L~P/»f M (x), 



X 



L 



" 7/ V x (x), 



(3a) 



(3b) 



where v, j3 and 7 are the usual critical exponents, and 
fi(x) are finite size scaling (FSS) functions with x = (T — 
Tc)L 1 / 1 ' being the scaling variable. Therefore, from the 
size dependence of M and x one can obtain the exponents 
/3/v and 7/V, respectively. 

The maximum value of susceptibility also scales as 
Jji/ V , Moreover, the value of temperature T* for which 
X has a maximum, is expected to scale with the system 
size as 



T*(L) 



(4) 



where the constant b is close to unity [23j. Therefore, the 
Eq. ((4]) may be used to determine the exponent 1/v. Wc 
have checked also if the calculated exponents satisfy the 
hyper-scaling hypothesis 



2/3/1/ + 7/1/ = D e 



(5) 



in order to get the effective dimensionality, D e ff, for both 
investigated AL lattices. 



III. RESULTS AND DISCUSSION 

The dependence of the magnetisation M on the tem- 
perature T, obtained from simulations on (3,4,6,4) and 
(3 4 , 6) AL with N = 6L 2 ranging from 384 to 98304 sites 
is presented in Fig. [21 The shape of magnetisation curve 
versus temperature, for a given value of N , suggests the 
presents of the second-order transition phase in the sys- 
tem. The phase transition occurs at the critical value Tc 
of temperature. 

In order to estimate the critical temperature Tc we 
calculate the fourth-order Binder cumulants given by Eq. 
(|2"c)) . It is well known that these quantities arc indepen- 
dent of the system size at Tc and should intercept there 

M. 
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FIG. 2: The magnetisation M as a function of the temperature T, for L = 8, 16, 32, 64, and 128 and for (3, 4, 6, 4) and (3 4 , 6) 
AL. 




FIG. 3: The susceptibility \ versus temperature T, for (3,4,6,4) and (3 4 ,6) AL. 



In Fig. |3] the corresponding behaviour of the suscepti- 
bility x is presented. 

In Fig. 4 the fourth-order Binder cumulant is shown 
as a function of the temperature for several values of L. 
Taking two largest lattices (for L = 64 and L = 128) we 
have T c = 1.590(3) and T c = 2.100(3) for (3,4,6,4) and 
(3 4 ,6) AL, respectively. 

In order to go further in our analysis we also com- 
puted the modulus of the magnetisation at the inflection 
M* = M(Tc). The estimated exponents j3/v values are 
0.180(20) and 0.103(7) for (3,4,6,4) and (3 4 ,6) AL, re- 
spectively. 

Basing on the dependence ln% on InL we estimated 
j/u = 1.46(8) and 7/1/ = 1.44(8) for (3, 4, 6, 4) and (3 4 , 6) 
AL, respectively. 

To obtain the critical exponent 1 /i>, we used the scaling 
relation (01 . The calculated values of the exponents 1/v 
are 0.83(5) for (3,4,6,4) and l/u = 0.94(5) for (3 4 ,6). 
Eq. ([3]) yields effective dimensionality of the systems 



D eS = 1.82(4) for (3, 4, 6, 4) and D eS = 1.64(5) for (3 4 , 6). 

The above results, indicate that the Ising 5 = 1 model 
on (3,4,6,4) and (3 ,6) AL does not fall in the same 
universality class as the square lattice Ising model, for 
which the critical exponents are known analytically i.e. 
P = i = 0.125, 7=| = 1-75 and v = 1. We have 
checked numerically, that Ising S = i model repro- 
duces these critical exponents with reasonable accuracy 
for both studied lattices [25[. We improved the value 
of the critical temperature Tc for these two lattices and 
S = \ as well, with respect to Ref. [2l| . 

The results are collected in Tab. fl] 

Except the exponent v, all critical exponents for 5 = 1 
differ for more than three numerically estimated uncer- 
tainties from those given analytically. 
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FIG. 4: The reduced Binder's fourth-order cumulant U as a function of the temperature T, for (3,4,6,4) and (3 4 ,6) AL. 



TABLE I: Critical points and critical points exponents for (3,4,6,4) and (3 4 ,6) AL. For comparison, the exact values for the 
square lattice Ising S = i model are included as well. 





s 


k B T c /J 


p/v 


7 A 1 


1/v 


D eB 


(3,4,6,4) 


1 


1.590(3) 


0.180(20) 


1.46(8) 


0.83(5) 


1.82(4) 


(3 4 ,6) 


1 


2.100(3) 


0.103(8) 


1.44(8) 


0.94(5) 


1.64(5) 


(3,4, 6,4) 


1 

2 


2.145(3) 


0.123(17) 


1.680(74) 


1.066(44) 


1.926(84) 


(3 4 ,6) 


1 

2 


2.784(3) 


0.113(10) 


1.726(8) 


1.25(13) 


1.952(22) 


square (4 4 ) 


1 

2 


2/arcsinh(l) 


i 

<s 


7 

4 


1 


2 
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